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O ■ Abstract 

We present a functional Schrodinger picture formalism of the (1+1)- 
dimensional 

• 0{N) nonlinear sigma model. The energy density has been calculated 

r'"! ' to two-loop order using the wave functional of a gaussian form, and from 

rS 

■ which the nonperturbative mass gap of the boson fields has been obtained. 

The functional Schrodinger picture approach combined with the variational 
technique is shown to describe the characteristics of the ground state of the 
nonlinear sigma model in a transparent way. 

ll.lO.Ef, ll.lO.Lm 

Typeset using REVTJt;X 



* e-mail : dkkim@phya . yonsei . ac . kr 
''^e-mail:ckkim@phya.yonsei.ac.kr 



1 



I. INTRODUCTION 



The ground state of the interacting quantum fields, in general, has a complicated struc- 
ture to investigate, rendering ordinary perturbation theories without much success. Espe- 
cially, concerning spontaneous symmetry breaking and bound states, perturbative ground 
states mislead to wrong results [T|]. In this respect, the functional Schrodinger picture (FSP) 
approach with variational approximation is expected to be a useful tool to examine the 
nonperturbative aspects of quantum field theory. 

In contrast to the usual perturbative expansion, the Schrodinger picture approach has 
the merit that one does not need to specify a particular Fock basis for the ground state of 
the Hamiltonian under consideration. Therefore, when there is no well defined Fock vacuum, 
this method appears to be a convenient choice 0. 

In particular, the nonlinear sigma (NLS) model [0-^ has a nontrivial vacuum structure 
that is composed of particle-antiparticle pairs. The ground state is not easily tractable 
applying the usual perturbation expansion. Therefore, it appears as a natural candidate for 
application of the FSP approach. The aim of this paper is to analyze the NLS model in the 
framework of the FSP method combined with a variational approach [^-Q, which is known 
to go beyond the perturbative scheme in some cases.. We will show that the nonperturbative 
phenomena like the mass gap and asymptotic freedom can be described in the Schrodinger 
picture in a direct way. 

The NLS model in lower dimensions has attracted much attention, since it has relevance 
to the low energy limit of QCD as well as condensed matter systems such as antiferromagnets. 
The NLS model in two dimensions is classically scale invariant and asymptotically free 
According to Mermin- Wagner and Coleman , the continuous symmetry can not be broken 
in 1 + 1 dimensions; the massless Goldstone bosons tend to acquire their masses. And it was 
shown that there is a mapping between the NLS model and the effective long wavelength 
action of the quantum Heisenberg antiferromagnet |[TT |. 

In section II, we briefly introduce the NLS model and its formulation in the FSP approach. 
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In the section III, we first calculate the energy density to two-loop order using a gaussian- 
type wave functional, and then derive the mass gap for the boson fields by minimizing the 
energy density. We will show that in the NLS model, the massive ground state is more stable 
than the massless one. In the last section, a brief summary and discussion of our results will 
be given. 



II. THE NONLINEAR SIGMA MODEL IN THE FUNCTIONAL 
SCHRODINGER PICTURE 



We start with the 0(iV)-invariant Lagrangian density [|,^ 



^ = ^d,<^ad^<^a, (1) 



where N scalar fields $a,a = 1, ■ ■ ■ , obey the constraint. 



N 



E = (2) 

a=l 

This constraint makes the theory quite complicated, since the N components of the scalar 
field $ are mutually dependent on each other. The coupling constant A is a measure of the 
strength of the self interaction of the scalar fields $a, and a small value of A corresponds 
to a weak interaction. The constraint in Eq.(2) means that one degree of freedom among 



the variables, is not a real dynamical variable |12]. Thus, we follow the standard 



prescription to get rid of the Nth field through the following nonlinear transformation 



- d, - - '^V^) 

(1 + 02/4)' '^ (1 + 0V4)' 
where 0^ = Y.a=i'Pa.'^- Substituting these expressions into the Lagrangian, we find the 
equivalent Lagrangian involving the — 1 fields (pa, 

^ 1 d^(f)ad^(l)a 

2A(l + </)V4)2- ^ 
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Here and in the following, the sum over the repeated indices are implied; otherwise, a 
comment will be given explicitly. 

The NLS model in terms of 0a's has no mass parameter, so these fields are classically 
massless Goldstone bosons. However, the Goldstone bosons originate from the breakdown 
of the continuous 0{N) symmetry, which can not occur in this case |T^. In the below, it 



will be seen that the Goldstone bosons become massive through the quantum mechanical 
self-interaction. 

The canonical quantization procedure of the classical system require the conjugate mo- 
mentum of the field (pa , which becomes 

dC 



7r„ 



1 



(5) 



A (1 + 02/4)2 

Thus, the Hamiltonian is written in terms of the canonical variables, tTq and (pa as 

where we used the conventions (p"^ = Y^a=i 4>a' ^ = Y,a=i ^"^^ 4>''^ = J2a=i {d-4'a/ dx)"^ ■ 

In the quantum theory, the dynamical field variable (pa{x) and TTa{x) become field oper- 
ators. These vTa's satisfy the canonical equal-time commutation relations with the 0a's such 
that 

[(pa{x), 7rb{y)]xo=yo = iSabS{x - y). (7) 

In the FSP representation of its quantum theory, the scalar field operator (pa{x) and its 
conjugate momentum T^a{,x) are realised as 



(Pa{x) (pa{x), 
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Now, we employ a Gaussian-type wave functional [§] with two variational function parame- 
ters Gab{x,y) and (j)a{x), 



1 



[det(27r^G)]i/4 



exp 



Gab i^^y) 



(My) - My)) 



(9) 



/x,y 4:h 

where repeated indices over a or b mean sums; otherwise, an explicit comment will be given 
in the below. This wave functional may be seen to obey that 

{^\M^)\^) = M^), 

(^Itt^x)!^) = 0. (10) 

Here, 0a(a;) is the expectation value of the field operator 4>aix), and the expectation value of 
the momentum operator has been chosen to be zero. From the above trial wave functional, 
we readily obtain the following results, which are needed for calculation of the expectation 
value of the Hamiltonian, 



{(f)aix)(pbiy)) = 0a(a;)06(y) + hGabix,y), 
{na{x)'Kb{,y)) = ^G-i^{x,y). 



(11) 



III. LARGE N CALCULATIONS 

In order to obtain the expectation value of the Hamiltonian with the gaussian trial wave 
functional it is necessary to handle the field operator in (1 + 0^/4)"^ properly. For 

that purpose, we expand it in terms of 0^/4, which involves an infinite number of terms of 
even powers of 4>a{x) fields. 

Thus, we have to evaluate terms as follows, 

(^|0^(x)^^~0^7r2(x)|^) and (^| 0^(x)~^^0^0'^(x)|^), (12) 

which can not be calculated in closed forms, except for some special limiting cases. Due 
to this inability, we require an approximation scheme. Here, we will resort to the large- iV 
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calculations [|T3| . The following example will show why the large- approximation is useful 
in the present problem; 



{'if\4>\x)4>\x)n\x)\'^) 



la4>a + hGaa){4>b4>b + hGbb)^Gj 
+ Tl 4>aGab4>bG^^ — -(f)a4>aSbb + —GabGbaG^^ — flGaa^bb 



(13) 



where the three terms in the brackets on the right-hand side are of the order of A^'^, iV^, and 
N respectively. Thus in the large limit, only the leading terms of A^^ order dominate. 



{^(t?{x)<j)\x)'K\x)\m) 



{(f)a(f)a + flGaa){Mb + hGbb)jGj, 

{^\(j)^{x)\^){^\(f)\x)\^){^\7r\x)\^). 



(14) 



Therefore it is clear that in the large A^ limit, the expectation values of the composite 
operators become |Q 

{^\(j)\x) ■ ■■(f)\x)n^{x){0T (j)'\x))\^) = {^\(f)\x)\^) ■ ■ ■ {^\(f)\x)\^){^\7T^{x){0T (j)'\x))\^) . 

(15) 



To construct an 1/A^ expansion [|T3| in a systematic way, we define a new parameter g such 
that 



g^XN, 



(16) 



where g is fixed to be finite. Thus, we are allowed to write the Hamiltonian expectation 
value in the following form 



2N 



1 + 



IN 
27 



1 + 



A/2 



(17) 



where on the right side, the expectation has been taken with respect to the gaussian wave 
functional in Eq.(9). Using the results for (0^), (vr^) in Eq.(ll), this equation can be rewritten 



as 
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I g 4 



(18) 



At this stage, we confine ourselves to the constant field configuration for the ipaix) fields, 
so that the square of the gradient of (j)a{x) in Eq.(18) vanishes. The scheme to expand 
(\l'|if |\E') in powers of h and discard terms higher than second order in will be adopted in 
the below. The % expansion is equivalent to the loop expansion |]5|,p!^. Thus, we are going to 
study the system to two-loop order. The energy density given in Eq.(18) can be expanded 
to second order in % to yield 



2 g 



2/(0) 

^ _ 1 hG, 



iGcc 



(V^^Cc), 



(19) 



2/(0) J 

where /(0) = 1 + 0^/4 has been defined. 

We now vary this equation with respect to Gbb{x, x) to determine the parameters Gbb{x, x) 
that minimize the energy expectation value. Thus, one gets the following relation: 



5{H) 



9 ^/(0) 



|G-i-/(0)[l + -A^a 1^ ' 
2 2/(0) ^bb 



G, ' 



5G, 



bb 



+ 



N h 

72/2(0) 



-[1 



2/(0) 



-GaaN' + 



n 



2/(0) 



SG 



bb 



0. 



(20) 



Here and in the below, we will use the following notations. The repeated indices over the 
letter b does not indicate a sum; while, the repeated ones over a means a sum. This equation 
gives the relation that the variational parameter Gbb must satisfy, 

4 3hGaa{z,z)\ _o 1 h 



Gbb (a;, y) 



1 



2/(0) 



6{x - y), 



21) 



9' n<P) V 2 /(0) 
which has also been expanded in powers of h and terms higher than h have been discarded, 
so that only terms up to can be retained in Eq.(19). Since it is practically impossible to 
solve this equation directly, we separate the equation into two parts as follows: 
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/ 



dp 



and 



TV 2 [l-{'ih/2f)Gaa{z,z)Y/^ J 27rVp' + m^ 



m = — — ^ — G^^[x,x), 



-,e-x.^[ip{x - y)] 



(22) 



(23) 



m 8 

where the latter one turns out to define the mass parameter of the boson operator (f){x) , as 
will be seen in Eq.(26). 

We analyze the equation using the following iterative method. First, we approximate 
the unknown function Gbb(x,y) by G^{x,y): 



g'Sm 



9 P{4>) f dp 



1 



=exp[ip(x - y)]. 



(24) 



N 2 J 2% Vp' + 

Second, to improve this approximation, we substitute the equation back into the coefficient 
of the right-hand side of Eq.(22). Thus, we have 



G^J{x,y) 



9 n<p) 



dp 



-^exp[ip{x - y)]. (25) 



N 2 [l-(3V2/)G'i°2(x,x)]V2 J 2n^p' + m'' 
Note that the multiplicative factor in front of the integral has a divergence involving a cutoff 
A. This must be removed by a proper renormalization of the wave function (j){x). Then, the 
wave function renormalized expression for Ghb{x,y) becomes 



G'd\x,y)^ 



9 P{4>) f dp 



I 



^exp[ip{x - y)]. 



(26) 



N 2 J 27r Vp' + 

This form will be used in the subsequent discussions in calculating the energy expectation 
value. As a result of this renormalization, Eq. (21) is now rewritten 



Gbb\^,y) 



-Vl + -r^G-^{z,z) 



S{x - y). 



(27) 



9^ f\4>)-^ ■ 2f($) 

Now, let us evaluate the Hamiltonian expectation value using Eq.(23) and (26). The 
V^Gaa in the Hamiltonian can be calculated multiplying Eq.(27) by Gbb{y, z) and integrating 
over the volume / dy: 



W^Gaa{x,x) 



9' r (0)^-1 



4 



Gaa {^^ ^) 



1 - ^^^Gaaix,x) 



2/(0) 



(28) 



The hamiltonian density has no classical potential energy part V{(f)), so we are allowed to 
set 

= 1- (29) 
Thus, the Hamiltonian to the second order in h, that is, to the two-loop order, is given by 

^j^^G-J{x,x)[l-^GUx,x)] (30) 

Using the mass defining relation Eq.(23) and two point Green's function Eq.(26), one 

can evaluate (if). 



N 



2w? 
9 


(1- 


^G 


la {X-, x) ) 


2m? 


[1- 


fig 


dp 1 


9 


8 J 


-A 27r ^/p^ + m 
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[1- 


fig 
167r 


M— )], 



(31) 

where A is an ultraviolet momentum cutoff. To make the energy density finite, we may 

renormalize the parameter of the theory. Defining the renormalized coupling constant gr, 

11 fi ^ ,4A2, 
g gr IGtt /i^ 

we obtain the finite energy density with the renormalization scale /x: 

(^^2^ ^ X). (33) 

Note that this is the result up to since the mass parameter is of fi order. This energy 
expectation value has a minimum away from the origin, since it is concave upward as 
increases. 

Hence, by minimizing the energy expectation value with respect to m^, one can derive 

the mass gap, 

d{H)/N _2_ n errf_ 
dm? gr Stt /x^ 

= 0. (34) 
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Prom this relation, one obtains the dynamically generated mass gap: 

IGtt" 



(m 



IJ, exp 



-1 - 



hgr 



(35) 



The NLS model has no dimensional parameters; the coupling g is dimensionless in two 
dimensions. But, we arrived at a dimensional parameter m^; this phenomena is an example 
of dimensional transmutation. At this value of the mass gap, the energy density to two-loop 
order becomes 



N 



9r 

h 



Stt 



exp 



-1 - 



1^' 
167r 

hgr 



(36) 



Note that the negative sign of the energy density indicates that the massive ground state is 
more stable than the massless one in the NLS model. 

We now return to the mass defining equation, Eq.(23). Differentiating it with respect to 
m^, one finds the relation between the coupling constant g and the cutoff A 

h f dp 1 



1 _h r 
g~8J 



27r ^/p^~+mF 



IGtt 



(37) 



This relation shows that as A becomes large, the coupling constant g approaches zero, thus 
satisfying the asymptotic freedom. This equation can be also rewritten in terms of the 
renormalized coupling gr and the renormalization mass scale in Eq.(32). 



IV. CONCLUSION 

In this paper, the functional Schrodinger picture approach has been applied to analyze the 
0{N) nonlinear sigma (NLS) model. We have considered the 0{N) NLS model in the large 
A^-limit and calculated the energy expectation value to the second order in /i(two-loop order) 
systematically, using the gaussian wave functional form. The Schrodinger picture approach 
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combined with the variational technique produced the mass gap and the asymptotic freedom 
of the ground state for the 0{N) NLS model in a straightforward manner. 

The majority of the literature on the 0{N) NLS model adopt the Lagrangian for- 
malism to investigate its nonperturbative phenomena, where a composite auxiliarly field 
cix) — J2a^a{x)^a{x) IS usually introduccd. However, here we discuss its nonperturbative 
phenomena directly without resorting to the superfluous auxiliarly field. 

The extension of our calculations to three dimensions will be straightforward.; the dif- 
ference is to write the Green's function in Eq.(26) in corresponding three dimensional form 
and carry out the integral in the three dimension. 
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